Optical properties and electronic structure of ZrBi2 
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We report optical (6 meV - 4 eV) properties of a boride superconductor ZrBi2 (Tc = 6 K) 
in the normal state from 20 to 300 K measured on high-quality single crystals by a combination 
of reflectivity and ellipsometry. The Drude plasma frequency and interband optical conductivity 
calculated by self-consistent full-potential LMTO method agree well with experimental data. The 
Eliashberg function a^^F^uj) extracted from optical spectra features two peaks at about 25 and 80 
meV, in agreement with specific heat data. The total coupling constant is \tr = 1-0 ± 0.35. The 
low energy peak presumably corresponds to the displacement mode of Zr inside B24, cages, while 
the second one involves largely boron atoms. In addition to the usual narrowing of the Drude peak 
with cooling down, we observe an unexpected removal of about 10 % of the Drude spectral weight 
which is partially transferred to the region of the lowest-energy interband transition (~ 1 eV). This 
effect may be caused by the delocalization of the metal ion from the center of the ,824 cluster. 

PACS numbers: 74.70.Ad, 78.20.Ci, 78.30.-j 



I. INTRODUCTION 

Boron, like carbon, covalently bonds to itself, easily 
forming rigid three-dimensional clusters and networks as 
well as planes, chains and even nanotubes. The discov- 
ery of superconductivity in graphite-like MgB2 at 40 K 
[H has stimulated intense research of other supercon- 
ducting boron phases, although they have so far shown 
quite modest transition temperatures. Zirconium dode- 
caboride ZrBi2 has one of the highest known Tc ~ 6 K 
among the binary borides except MgB2. Even though 
the superconductivity there has been discovered almost 
40 years ago , only recent progress in the single-crystal 
growth enabled its extensive studies. 

In dodecaborides MB12, boron atoms make up a three- 
dimensional network forming spacious B24 cages which 
accommodate metal ions. The isotope effect in ZrBi2 for 
zirconium i (/? « -0.32) is much larger than for boron 
[3| (/? ~ -0.09), pointing to a large contribution of lat- 
tice modes involving Zr atoms to the electron-phonon 
coupling responsible of superconductivity in ZrBi2. The 
inversion of specific heat measurements gives a pro- 
nounced peak in the phonon density of states at about 
15 meV, attributed to Zr vibrations in oversized boron 
cages. The coupling to this mode can also be enhanced 
by strong anharmonicity. The contribution to electron 
phonon coupling of this mode was also reported in See- 
beck effect measurements @. 

Notably, significant controversy exists regarding the 
strength of the electron-phonon interaction. According 
to the McMillan formula, a weak-to-medium coupling 
A ~ 0.68 provides the observed Tc, given the phonon 
mode frequency of 15 meV and a 'standard' value of 
the screened Coulomb potential ^jl* ~ 0.1 — 0.15 Q. 
A weak coupling was confirmed by specific heat mea- 



surements @. However, tunneling Q and point-contact 
spectroscopy yield 2Ao/fcs'rc « 4.8 suggesting a very 
strong coupling regime. This discrepancy between the 
results of surface and bulk probes has been ascribed to 
the surface-enhanced superconductivity in ZrBi2 [1, [§]. 
In principle, the value of A can be also derived from the 
DC resistivity measurements, although this calculation 
relies on the unknown value of the plasma frequency. 

Band structure calculation fio'l shows that B 2p and 
Zr 4d states contribute almost equally to the density of 
states at the Fermi level. This makes ZrBi2 electroni- 
cally very different from MgB2, where Mg states have a 
vanishing contribution to the metallic and superconduct- 
ing properties pd| . However, no experiments which can 
directly verify the prediction of the band theory, such 
as angle-resolved photoemission (ARPES), de Haas- van 
Alphen (dHvA) effect, have been done so far on this com- 
pound. 

Optical spectroscopy that probes free carrier charge 
dynamics and provides experimental access to the plasma 
frequency, electron scattering and interband transitions, 
may help to clarify a number of existing uncertainties 
and test the predictions of the band calculations. By 
a combination of reflectivity and ellipsometry, we have 
obtained spectra of optical conductivity and dielectric 
function in the broad range of photon energies and fre- 
quencies. We concentrate on general optical fingerprints 
of the electronic structure of ZrBi2 in the normal state. 

For comparison, we performed calculations of the band 
structure and corresponding optical properties. The 
electron-phonon coupling function a1^F{ui) has been ex- 
tracted from the fit of optical data and compared with 
thermodynamic and Raman measurements. 

An unusual decrease of the Drude plasma frequency 
with cooling down is observed. By comparing with other 
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FIG. 1: (color on-line) Crystal structure of ZrBi2. The clus- 
ters Bi2 and B24 are displayed as polyhedrons. The unit cell 
containing 13 atoms is colored in red. 



clustered compounds, we found that the delocaUzation 
of the metal ion (weakly bonded to the boron network) 
from the center of the boron cage could be at the origin 
of large changes in the electronic properties. 



II. CRYSTAL STRUCTURE 

ZrBi2 crystallizes in the UB12 type structure jT3| which 
can be viewed as a cubic rocksalt arrangement of Zr and 
B12 cuboctahedral clusters (Fig[T|). For our discussion, 
it can be more convenient to represent this structure as 
a face centered cubic structure of Zr, surrounded by B24 
cages. 



III. OPTICAL EXPERIMENT AND RESULTS 

Large single crystals of ZrBi2 were grown using a high 
frequency induction zone furnace [l3l |. Optical measure- 
ments were performed on (001) surface with dimensions 
4x4 mm^ . The surface of this extremely hard material 
was polished using diamond abrasive disks with finest 
grain size of 0.1 /xm. In the photon energy range 0.8-4 
eV, the complex dielectric function e{cj) was determined 
directly using spectroscopic ellipsometry at an incident 
angle of 60°. For photon energies between 6 meV and 
0.8 eV, the reflectivity of the sample was measured using 
a Bruker 113 Fourier transform infrared spectrometer. 
For both ellipsometry and reflectivity experiments, the 
sample was mounted in a helium flow cryostat allowing 
measurements from room temperature down to 7 K. The 
reference was taken by in situ gold evaporation. 

Fig [2^ presents the real and imaginary parts of e{Lu) 
measured by ellipsometry at selected temperatures. The 
reflectivity measured at low frequencies and extracted 
from the dielectric constant in the visible range, are plot- 
ted in FigHs. 

In order to obtain the optical conductivity in the in- 
frared region we used a variational routine [lj| yielding 
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FIG. 2: (color on-line) a) Real and imaginary part of the 
dielectric function measured by ellipsometry. b) Reflectivity 
measured in the infrared and calculated from the dielectric 
constant at higher frequencies. Inset: low frequency region. 



the Kramers-Kronig consistent dielectric function which 
reproduces all the fine details of the infrared reflectivity 
data while simultaneously fitting to the complex dielec- 
tric function in the visible and UV-range. This procedure 
anchors the phase of the infrared reflectivity to the phase 
at high energies measured with ellipsometry [TH - 

Fig. [3^ shows the evolution of the optical conductivity 
with temperature. One can see a Drude-like peak that 
indicates a metallic behavior. The DC conductivity (ctoc) 
@, 01 plotted as symbols in the inset of Fig[3^, agrees 
reasonably well with the extrapolation of the optical data 
to zero frequency. 

FigHb depicts the effective number of carriers 



2mVr 
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where m is the free electron mass, e is the electron charge, 
Vc=101.5 A"^ is the unit cell volume. We observe a trans- 
fer of spectral weight to low energy, which comes from 
the narrowing of the Drude peak. For higher energies 
(above 100 meV) an unusual decrease of up to 10% is 
observed. This will be discussed in section IV CI 
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FIG. 3: (color on-line) a) Optical conductivity of ZrBi2 at 
selected temperatures. Inset: low energy part, symbols repre- 
sent anc from 'i1].The dotted line shows the interband conduc- 
tivity from LDA calculation, b) Effective number of carriers 
per unit cell of 13 atoms, inset: low frequency part. The 
horizontal dotted line corresponds to the calculated plasma 
frequency of 6.3 eV. 



IV. LDA CALCULATIONS 
A. Band structure 

The band structure was calculated for the unit cell con- 
taining 13 atoms. We used a full potential plane waves 
Linear Muffin-Tin Orbital (LMTO) program ^ within 
a Local Spin Density Approximation (LSDA) and Gen- 
eralized Gradient Approximation (GGA) [T?!. The band 
structure for high symmetry directions in the Brillouin 
Zone (BZ) is shown in FigH^. The density of states 
(DOS) calculated by the tetrahedron method is displayed 
in Fig[3)D. It appears clearly that the main contribution 
to the DOS at the Fermi level comes from zirconium 4d 
and boron 2p states. The total DOS at the Fermi level 
N{Ef) = 1.59 eV"icell~^ N{Ef) for zirconium d elec- 
tron is 0.55 eV~"'^cell~^ and the one for p electrons of B12 
clusters is 0.65 eV^^cell^^ which is in good agreement 
with a previous report [10]. 

Mapping the repartition of charges in a crystal is help- 
ful to understand the chemical bonding. The calculated 
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FIG. 4: (color on-line) a) Calculated band structure, b) par- 
tial density of states (PDOS) of ZrBi2. c) Charge density of 
in the a-b plane of the cubic cell. 



charge density of valence electrons in the a-b plane of 
the cubic cell is displayed in Fig [It. The high density 
between direct neighbor boron atoms indicates a strong 
covalent bonding while the ion in the middle of the B24 
cage (Zr in our case) is weakly bonded to the boron skele- 
ton. This has strong consequences on vibration modes in 
the crystal and may explain the presence of a very low 
frequency phonon (~ 15 meV) associated with a displace- 
ment of the zirconium atom in the B24 cage. 

FigOi shows the Fermi surface in the first Brillouin 
zone. The sheet centered at X (marked in red in Figl5]) 
has electron character while the network centered on the 
r point (marked in blue in FiglS]) is of hole character. 
The two sheets are touching at point a along the F X 
direction. In contrast to MgB2, which exhibits a strongly 
anisotropic Fermi surface composed of different sheets 
with very different character of carriers, the Fermi surface 
in ZrBi2 is much more isotropic. 



B. Optical properties 

The LDA band structure allows us to calculate the 
optical parameters. The interband contribution to the 
real part of the optical conductivity can be deduced from 
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FIG. 5: (color on-line) Fermi surface of ZrBi2 containing two 
sheets: a hole-like network (blue) and electron-like closed sur- 
faces (red). In the left figure only the first sheet is displayed 
for clarity. The points of tangential contact between the two 
sheets are designated by symbol a. 



the equation: 



CTi (w) 



V / d^k\Pf,\^5[Ef{k)-E,{k)-nw] 



Here, Ei{k) and Ef{k) are the energies of the initial (oc- 
cupied) and final (empty) states, respectively, k is the 
wave vector inside the BZ where the transition Ei{k) — > 
Ef{k) occurs. Only direct transitions are taken into 
account. The calculated optical conductivity is shown 
as the dotted line in Figl3^ together with experimental 
curves. At frequencies above 2 eV a very good agree- 
ment is observed. At lower frequencies, the experimental 
curves are systematically higher likely due to the tail of 
the Drude peak. 

The plasma frequency was computed through integra- 
tion on the Fermi surface : 



^^p,LDA 



127r2 
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VfdS 



(2) 



where Vf is the Fermi velocity. We found a value 
^p.LDA = 6.3 eV which corresponds to N^fj ~ 2.9 (dot- 
ted line in Fig[3lD). The experimental value of Neff{Lu) 
reaches this level at about 2 eV where the experimental 
optical conductivity (Fig[3^) shows a minimum. This in- 
dicates that the calculated plasma frequency is close to 
the experimental one. 



V. ANALYSIS AND DISCUSSION 
A. Extended drude analysis 

Fig [6] shows the frequency dependence of the mass 
renormalization m* (w) / mi, {mi, is the band mass) and the 
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FIG. 6: (color on-line) Extended Drude analysis of the optical 
conductivity of ZrBi2 for different temperatures. Inset shows 
the T dependence of the plasma frequency used to calculate 
m*/m ans I/r. 



scattering rate 1/t(cj) obtained by the extended Drude 
formalism: 



m*{uj) 
mh 
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The plasma frequency fip was obtained at each temper- 
ature by intergation of optical conductivity from zero to 
0.74 eV (6000 cm^^) (inset in FigO. This cutoff fre- 
quency was chosen such that fip is close to the plasma 
frequency obtained using the fitting procedure (fit 3) de- 
scribed below. The temperature dependence of fJp is too 
large to be neglected as it is usually done. 

At low frequencies, the decrease of the scattering rate 
with cooling down indicates a narrowing of the Drude 
peak. The strong frequency dependence below 100 meV 
of 1/T(a;) and m*(a;)/mf, indicates a significant electron- 
boson, presumably electron-phonon interaction. 
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B. Electron-phonon coupling 

The signature of electron-phonon interaction shown by 
the extended Drude analysis makes the simple Drude 
model inapplicable to describe the low frequency region. 
Therefore, we adopted the following model for the dielec- 
tric function: 



uj[io + iM{Lu,T)] 



E 



^2 . 



(5) 

Eoo represents the contribution of core electrons, the sec- 
ond and the third terms describe free carriers and in- 
terband contributions respectively. The latter is taken 
to be a sum of Lorentzians with adjustable parameters. 
The frequency dependent scattering of the free carriers 
is expressed via the memory function: 



M{uj, T) = 7,™p - 2i / dnalF{n)K 



ui 



2ttT ' 2nT 



(6) 

where jimp is the impurity scattering rate and af^Fj^}) is 
the transport Eliashberg function. The kernel K is [I^ : 

Kix,y) ^ 1 + {^[^-(1 -ix + ly) - *(1 + ly)]] - 

{y ^ -y} (7) 

where ^(x) is the digamma function. 
The reflectivity at normal incidence is: 



R{lo) 



(8) 



We fit simultaneously R{lo) at low frequencies and both 
ei{uj) and 62(0;) at high frequencies. 

As a starting point, we used a'^F{il) deduced from 
specific heat and resistivity measurements (neglecting 
the difference between isotropic a'^F{Vt) and transport 
af^Fip)) and adjusted all other parameters (fit 1: 
solid red curve in Fig [7^). This function (solid curve in 
FiglH]) has two peaks around 12 meV and 60 meV and 
gives A = 0.68. A clear improvement of the fit quality 
was achieved compared to one with a single Drude 
peak (dashed dotted curve) with frequency independent 
scattering rate. The refiectivity curve is within absolute 
reflectivity error bars of 0.5%. This means that optical 
data are consistent with the mentioned value of the 
coupling constant and the observed Tc [1, Q ■ 

It is also important to determine the spread of param- 
eters consistent with optical data alone. It is known that 
the extraction of a'^F{Vl) from optical data is an ill-posed 
problem [H, 113, El)] which means that fine details cannot 
be extracted. We model the Eliashberg function as a 
superposition of Dirac peaks: 



E 

k 



Ak5{n-xk) 



(9) 



TABLE I: Set of parameters used in fit 2 and fit 3. In fit 3, 
the Lorentzian L\ is split into three difi'erent ones. 
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and treat Ak and Xk as free parameters. We did not 
succeed to obtain a reasonable fit with a single mode 
even using some gaussian broadening within 10 meV but 
two peaks are sufficient to reproduce the shape of exper- 
imental curves. The central frequencies where found to 
be a;i=208 ± 19 cm'^ and X2=602 ± 128 cm^^. The 
weights are ^i=83 ± 21 cm"^ and ^2=94 ± 35 cm~^ 
respectively. The electron-phonon coupling constant can 
be derived from the parameters of the fit by the relation: 



Xtr — 2 



dna^^F{Q) 



= 2y^ 



(10) 



In our case, this ratio leads to a partial coupling constant 
At,.,i=0.82 ± 0.28 for the mode at xi and Atr,2=0.18 ± 
0.10 for the mode at X2- The total Atr=1.00 ± 0.35. The 
other parameters of the fit are displayed in table I. The 
reflectivity corresponding to this fit (fit 2) is displayed as 
the dotted (green) curve in Figl?^. 

The low frequency mode can be logically assigned to 
the vibration of weakly bounded zirconium ion inside B24 
cages, although the ab-initio calculation of the phonon 
structure is needed to substantiate this assignment. The 
high frequency mode presumably corresponds to a rigid 
boron network vibration; a situation similar to interca- 
lated fuUerenes and other clustered covalent structures. 

The error bars of Xk and Ak are quite large because 
of a significant numerical fit correlation between the low- 
est energy interband transition and the parameters of 
the oscillator (LI in table I) and those of af^F{n). In 
other words, it is difficult to decouple the effects of low 
energy interband transition and electron-phonon inter- 
action on the optical conductivity. Therefore we per- 
formed another fit (fit 3) where the interband contri- 
bution was modeled (parameters in table I) to match 
the one of the interband optical conductivity calculated 
by LDA (solid line in FiglTja). The presence of sharp 
structures slightly worsens the quality of the fit. The 
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FIG. 7: a) DifTerent fits of reflectivity at 25K. dasli-dotted 
line is a fit witfi a single Drude peak. Fit 1 with a'^F{Q,) from 
01 (solid red curve). Fit 2 and 3 atrF(0.) was adjustable 
as described in the text. b)Interband contribution to optical 
conductivity calculated from LDA (solid curve) and obtained 
from the Fit 3 as described in the text. 



low energy mode of a^^F{fl) is almost unchanged by this 
procedure (xi=227 cm~^, 74i = 106 cm~^ and Atri=0.93) 
but the second mode tends to shift to higher frequencies 
(x2 = 1003 cm-\ v42=190 cm~i and Atr2=0.38) and the 
total coupling constant becomes Atr=l-3. The reflectiv- 
ity corresponding to fit 3 is plotted as a dashed blue curve 
in Fig[7^. 

The central frequencies Xk and partial electron-phonon 
constants are shown as symbols in Figl5]for fit 2 (cir- 
cles) and fit 3 (squares) described above. The Eliashberg 
function from Q used in fit 1 is plotted as a solid curve. 
Central frequencies of the phonons are larger than those 
derived from specific heat and resistivity measurements 
[| (140 cm-i and 420 cm^^) (Fig. ^. 

A more traditional way to extract ai^i^(ri) from 
optical conductivity is to use directly the inversion 
formulafillii: 
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FIG. 8: (color online) The solid curve (fit 1) corresponds to 
a'^F{Q,) from i5|. Central phonon frequencies and partial cou- 
pling constant extracted from flt of optical data with fit 2 (ci- 
cles) and flt 3 (squares) procedures as described in the text. 



One can also calculate the integrated coupling constant 



(12) 



(11) 



Taking the frequency derivative twice requires heavy 
spectral smoothening of optical data. The result of inver- 
sion is given in Figl9]for different frequency smoothening 
windows. The robust features in af^F{u>) are the peaks 
at about 250 cm~^ and 900 cm~^, in good agreement 
with the results obtained by fitting procedures described 
above. Since Eq. fTTjl assumes the absence of low en- 
ergy interband transition, in our case some spurious fea- 
tures in a'^j,F{u!) may appear at frequencies above ^^1000 
cm~^. Even though smoothening has a significant in- 
fluence on the shape of afj.F{uj), all curves Xtr,int{^^) fie 
almost on top of each other and give a value of about 
1.1 around 1000 cm~^. Summarizing the results of all 
previously described approaches to extract a^^F^to), we 
obtained Xtr = 1-0 ± 0.3. This value indicates a medium 
to strong coupling in ZrBi2. 



C. Temperature dependence of Drude spectral 
weight 

The plasma frequency flp at 25 K, obtained from differ- 
ent flts (Table I) is 5.2-5.5 eV. The uncertainty is coming 
from the fact that interband contributions, which can not 
be separated unambiguously from the Drude peak. This 
value is slightly lower than the one calculated by LDA 
(6.3 eV). However this agreement is good and it sug- 
gests that LDA describes well the conduction bands near 
the Fermi surface. The temperature dependence of the 
plasma frequency though appears to be anomalous. As it 
was mentioned above, the integrated number of carriers 
decreases for energies above 100 meV when temperature 
is decreased. This is shown in FigfTUl where we plot N^ff 
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FIG. 9: (color online) ^^^^(fi) obtained using Eg. ([11) at 
T=25 K (solid curves) and integrated coupling constant 
Mr,int(io) (dashed curves) for different smoothening windows. 

as a function of temperature normalized by its value at 20 
K for selected cutoff frequencies. The effect is maximal 
(about 10%) at around 0.7 eV. Normally, the narrow- 
ing of the Drude peak with a temperature independent 
plasma frequency results in a slight (1-2%) increase of 
Neff at frequencies 0.5-1 eV. In our case one has to as- 
sume that the plasma frequency decreases with cooling 
down. On the other hand, spectral weight is almost re- 
covered at around 2 eV which is due to the fact that it 
is transferred from the low frequency range to the region 
between 0.5 and 2 eV (FigO]). 
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FIG. 10: (Color online) Temperature dependence of relative 
effective number of carriers for different cutoff energies. 

We do not have a definitive explanation of this effect. 
In order to check the direct effect of temperature on the 
band structure, we have performed LDA calculations at 
finite temperatures by adding a Fermi-Dirac distribution 
for the bands filling in the self consistent calculation. 
The thermal expansion was also taken into account using 
data from Ref,5;. We found that the calculated plasma 
frequency decreases around 0.4% with increasing tem- 
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FIG. 11: a) Representation of displacement directions in the 
-624 cage, b) Plasma frequency versus ionic displacement 
along [100], [110] and [111] directions. 

perature from K to 600 K. Both sign and amplitude of 
the change in the calculated plasma frequency cannot ex- 
plain the experimental observations. The spectral weight 
due to interband transitions has also been calculated and 
found to be almost temperature independent. Therefore 
a decrease of 10% of the spectral weight is so big that it 
likely results from a significant change of the electronic 
band structure. 

One can speculate that a certain structural instability 
is due to the ability of zirconium atoms to move almost 
freely in spacious boron cages. According to Refl23l. a 
displacement of the metal ion in the i?24 cage of about 
5% of the unit cell is responsible for the symmetry break- 
ing which makes low frequency Zr mode Raman-active. 
We have therefore calculated the band structure and the 
plasma frequency as a function of displacement of the 
metal ion in the directions [100], [110] and [111]. The de- 
localization configurations are sketched in FigfTTk. 

Electron states at the Fermi surface are equally coming 
from Zr atoms and B12 cluster (FigHjD). Therefore the 
highest plasma frequency is obtained when Zr ion is in the 
center of the B24 cluster. The plasma frequency decrease 
will be more pronounced when Zr shifts towards a region 
with lower charge density (Figlllb). 

Such an anomalously large atomic displacement has 
already been observed in clathrate-type systems that 
present a similar cluster structure [U, [2^ • In these com- 
pounds, the displacement of an ion located in the middle 
of a cage is strongly anisotropic and can be described by 
a fractionally occupied split site. It was also found that 
this atomic disorder was present from room temperature 
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down to lOK. 

In the picture of a multiple well potential, the delo- 
calization due to thermal excitation tends to center the 
ion in the middle of the cage. At low temperature, the 
ion is located in the bottom of the well out of the central 
position. Further structural investigations are necessary 
to elucidate the relevancy of this or other scenarios in 
dodecaborides. This kind of spectral weight transfer is 
also predicted by the model where the charge carriers 
are small Dolaronsp6|. While keeping in mind that they 
have been derived for a low charge carrier density, the 
observed direction of the spectral weight transfer agrees 
with those models. 



VI. CONCLUSIONS 

In this work, optical properties of zirconium dodeca- 
boride have been studied as a function of temperature. 
We found that LDA theory predicts correctly interband 
optical conductivity and free electron plasma frequency. 

ZrBi2 appears to be a good metal with a bare plasma 
frequency of about 5.5 eV. The Eliashberg function 
a'^j.F{fl) was extrated by two methods giving consistent 
results. Two peaks where found in aj^F(ri) at about 
200 and 1000 cm~^ which is consistent with specific heat 
measurements |5||. A medium to high electron-phonon 
coupling regime is found with A^^ ~ 1.0 in agreement with 
a high ratio 2Ao/kTc = 4.8 from tunneling measurements 



. The low frequency peak that likely corresponds to a 
vibration mode of zirconium in boron cages is the main 
contributor to the coupling constant. 

However, this compound shows deviations from con- 
ventional metallic behavior. Namely, a significant 
anomaly in the temperature evolution of the integrated 
spectral weight was observed. About 10% of the spectral 
weight is removed from the Drude peak and transferred 
to the region of low energy interband transitions. The 
origin of this anomaly has not been clearly identified. 
However, our simulation of the impact of the delocaliza- 
tion of the metal ion in the boron cages shows that it 
could be responsible for the large change in the plasma 
frequency. 
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